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ABSTRACT 


The turning characteristics of the Bell 100 ton captured 
Air Bubble Surface Effect Ship (CAB/SES) are studied with 
the aid of the Oceanic Incorporated computer simulation. 
The craft's motions for straight runs, straight runs with 
waves and turning runs under calm water conditions, are 
compared with the final set of computer simulation turning 
maneuvers with waves involving changing speeds, wavelengths, 


wave amplitudes, and final rudder angles, 
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thrust, also tangential force 

longitudinal velocity 

lateral velocity 

vertical velocity 

horizontal distance in direction of forward motion 
longitudinal force 

horizontal distance directed to starboard 

lateral force 


vertical distance directed downward, also 
heave displacement 


vertical distance from keel to center of gravity 
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mass density 
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I. INTRODUCTION 


- 


А. CAB/SES DEVELOPMENT 

The very latest design in displacement ships still 
manage to operate very inefficiently in the 40-knot-plus 
range. In fact the state-of-the-art in hull design peaked 
in the early part of this century. Ship designers have 
realized all along that ordinary displacement vessels 
sacrificed much efficiency for added thrust at speeds in 
the 40 knot range due to hydrodynamic drag. Propulsive 
thrust required increases approximately as the square of 
the desired velocity. To bridge this 40 knot gap, naval 
architects turned toward the surface effect ships (SES). 
Hovercraft, hydrofoil, air cushion vehicles (ACV), ground 
effects machines (GEM), and captured air bubble (CAB) craft 
all come under this generic grouping. A SES is basically 
any ocean going or amphibian craft which has the capability 
to operate in the air/water interface, thereby greatly 
reducing hydrodynamic drag. This study utilizes the CAB; 
in particular, the BELL 100B CAB/SES. 

The CAB was chosen over other SES due to projected 
state-of-the-art feasibility for high speed and sea keeping 
capabilities. Whereas hydrofoils and ACV's are size limited, 
theeCAB Ys net. The CAB is mainly supported by an air 
bubble enclosed in the hard structure of a deck, two 
sidewalls, port and starboard, and flexible seals at the 


bow and stern. The air bubble is provided by a blower fan 
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system which provides enough pressure to support and lift 

the craft. Unlike the ACV cousin, the CAB retains its air 
pressure with the exception of leakage and occasional 

venting in higher sea states and maneuvers. The craft 

reacts much like a displacement craft when not supported on 
the bubble but maneuvering is extremely sluggish. Propulsion 
in several CAB/SES already built has been either waterjets 

or supercavitating propellers. 

When "on-the-bubble", there is a pocket or depression 
in the water equivalent to the gross displacement minus that 
part of the displacement which is being supported by the 
Sidewalls and seals. As the craft increases speed and moves 
forward so too does the depression. This causes the same 
type of waves as a displacement hull. AS acceleration 
continues, the CAB overrides the forward portion of the 
depression. This is called the secondary hump and a bow-up 
trim is obtained. Continuing acceleration causes the 
depression to slip aft under the stern seal so that a smooth 
flow results. Hump speed has been surpassed, wave drag and 
skin friction drag have been significantly diminished, and 
the drag E velocity curve decreases to values well 
below drag experienced by displacement vessels. 

Over the past few oem a joint research program has 
been carried out by the Navy and the Maritime Administration 
under the heading of Joint Surface Effects Ship Project 
Office (JSESPO). Presently that office is known as SESPO 


Since the project is now sponsored by the Navy alone. 


dad: 





Several craft have since been designed, built, and tested. 
Among these are the XR-3, which is presently undergoing 
further research and evaluation at the Naval Postgraduate 
School in Monterey, California. Two 100 ton craft, the 
Aerojet General 100-A, propelled by a gas turbine water jet, 
and the Bell 100-B, a supercavitating propeller driven craft, 
are also undergoing evaluation by the SESPO. The ultimate 
projection of this program is to provide a CAB/SES in the 
multithousand ton class with cruise speeds from 80 to 100 
knots. 

The next step the project is expected to take is con- 
struction of a CAB/SES in the 2000 ton range. However, 
there are presently many questions such as seakeeping 
ability, safety, and dynamic behavior, which must be 
answered. Of the several methods for answering questions 
such as these, the most expensive would be to construct the 
full size prototype. On the other hand, towing tank model 
trials could be run, but probably the most effective method 
would be a computer simulation. In the long fun the computer 
Nacion data would have to be corroborated by empirical 
data. If the simulation were exceptionally good or at 
least a close approximation perhaps millions of dollars and 
many manhours could be saved. In the case of the CAB/SES 
the state-of-the-art rests in the two 100 ton craft already 
built. Instead of immediately building a 2000 ton craft 
based upon minimal data already obtained, Oceanics Incorpo- 


rated was contracted to provide a computer simulation, for 


О 
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the loads and motion study, see computer program section for 
the Bell 100B, with the following overall objectives: 


J. To provide a computerized form of the six degree 
of freedom equations for SES craft that will 
yield time domain outputs of motions for the on- 
bubble mode of operation. 

2. To determine the motions of specified SES confi- 
gurations with special reference to human habit- 
ability, platform requirements (for military or 
commercial applications) and maneuvering require- 
ments for the on-bubble mode of operation. 

3. To establish a mathematical model for predicting 
maximum wave induced bending moments, both 
longitudinal and transverse, on SES craft to 
provide techniques that would yield slamming 
response for the on-bubble mode of operation. 

Ц. То determine the stability derivatives in the 
vertical plane from measured data of the 
vertical plane motions of existing SES craft 
in a random seaway. 


Direct comparison of simulated data with experimental data 
should then provide a reasonable measure of the accuracy 
of the simulation. The next logical step is to convert 


the 100B program to simulate a 2000 ton craft. 


B. BELL 100B CAB/SES PROGRAM 
J P Ol mulation Statistics 
"All statistics given were held constant throughout 
ten but may be modified. Parameters listed as 
variable were in fact varied throughout the procedure. 
Total Craft Weight: 209,999.3 lbs 
Center of Gravity: 


kengpstcndinal (fe fwd. Transom): 33.09 ft 
Vertical (ft above keel): 7.09 ft 





я Kaplan, Paul, James Bentson, and Theodore P. Sargent, 
Scud ok Surface Effect Ship (SES Moti 
[am К ations of Motion of SES Craft wrth Six Degrees of 
reedom, p. 2, prepared for Department o e Navy Suriace 
Effect Ships Project Office, Contract No. C-1-35055, August 1971. 
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Total Craft length: 72 ft 


Fetal Craft width: ВО ВЕ Кросс сана starboard 
sidewalls) 


Ventral Fins: ВЕ У 3 16 
Plenum: | 
Length waterline: 65.31 ft 
Length deck: 65.31 ft : 
Width waterline: BIOS 
Width deck: 28.00 ft 
Height average: Celie tt 
Propulsion: Twin screws variable thrust 
Steering: Twin rudders 
` Turn Rate: Variable 
Maximum: Variable 
Speeds: Variable 
Fans: 


Bowseal: 2 
Sternseal: 2 


о. MDoscriptron or FErogram 

The Bell 100B computer simulation is a Fortran IV 
program primarily designed to calculate time histories of 
motions and loads of a surface effect ship, operating in 
the on-the-bubble mode in a seaway in six degrees of freedom. 
The program is of modular form in that each major element 
of the craft has its own subprogram, 18 in all. This feature 
enables the programmer to make a wide range of parametric 
changes by simply changing input data. Additionally, this 
modular format facilitates changes or additions to any 
part of the program. 

For any given program run, the en program calls 
the input and initialization subroutine INCON. The program 


then calculates the motion time histories by calling the 
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various Modular subroutines for as many time increments as 
specified to complete the test case. During the calcula- 
tions, the subroutines will, 1f requested, output variables 
selected by print option switches at desired print intervals. 

Reference (1), the users manual, contains a break- 
down on the contents of each subroutine. However, since 
the all important input data is handled by INCON, INCON will 
be discussed.  INCON is composed of 19 separate blocks. 
Within each block there are one or more options or option 
control tags. Note in oE I, some of the options are 
listed as not required. The other options, those required, 
are considered to be the basic or specific data options 
without which the program will not run. 

The original program's output occurred on a print 
time basis dependent on which print option switches were 
utilized. Additionally in the original program, the user 
could request print options which produced vertical and 
lateral plane summaries. See the computer program below 
as modified for this study. 

3. Formulation of State Equations of Motion 

The equations of motion utilized by the simulation 
are formulated in terms of the various force or moment 
contributing elements. 


n ; 
F = у Ер (1) 


Total iel 1 


Jus 





Table I. CAB/SES Incon Data Options 


BLOCK OPTION REQUIRED DESCRIPTION 


NUMBER 
001 01 Yes Timing information 
02 No Print options 
03 Yes Integrators and error tolerances 
04 No Summary prints 
05 Мо Motion option switches 
002 01 Yes* Summary mass properties 
02 Yes* Discrete mass properties 
003 00 Yes Craft geometric input 
004 01 Yes Sidewall input 
02 Yes Appendage input 
005 00 Yes Sternseal input 
006 00 Yes Bowseal input 
007 - 01. . . Yes Plenum geometry 
02 Yes Froude number input 
008 01 Yes Propulsion input 
02 No Engine out input 
009 01 Yes Rudder input 
02 No Rudder motion input 
010 00 Yes Aerodynamic input 
011 01 Мо Waves; amplitude vs. frequency 
02 No Waves; amplitude vs. wave length 
012 00 Yes Initial condition 
013 00 Yes End of case indicator 
014 00 Yes End of run indicator 
015 01 No Shears and moments 
02 No Shears and moments 
03 No Shears and moments 
016 Not Used 
017 Not Used 
-018 00 No Case title card 
019 01 Yes Bowseal fans 
02 Yes Cushions fans 
03 Yes Sternseal fans 


*either but not both of these options 


16 





Fi " "bubble, bi “sidewalls, 7 "ρου seal, S “rudder, 


+ F + F 
waves 


propulsion, š e rodante 5 


Í 


* Fstern seal, E "gravitational, "cs 
where F represents a moment or a force and 1 represents a 
| X sa ре such as in the lateral direction or 
| round the Larsen. PO ud adi poe 

The craft's trajectory is referred to an inertial 
reference frame, a right handed cartesian coordinate system. 
Thus, Хо, У; апа 2 show the trajectory with the craft's 
center of gravity at time t = 0 as the point (0,0,0). Roll 
($), pitch (8), and yaw (y) angles are also referenced to 
this inertial coordinate system. Displacements toward the 
bow are considered as in the positive Xo direction; displace- 
ments to starboard are considered as positive in the y 
direction; ο... vertically down from the center of 
gravity are considered to be in the positive 2 direction; 
yaw (Y) is positive with a turn to starboard; roll is 
positive with a heeling motion to starboard; and a bow up 
attitude is a positive pitch angle. Normal ship's motions 
including the CAB/SES have small values of pitch (8) and 
roll (¢). The main angle indicating the orientation of 
the craft relative to the inertial frame is yaw angle (y) 


which may have a large range if maneuveuring occurs. Motion 


ДЕ 





in the 2 direction, heave, will also be greatly limited. 
Since roll, pitch and heave values remain relatively limited 
in magnitude, overall complexity in simulating hydrostatic 
and hydrodynamic forces and moments jis greatly reduced. 

The axis system describing the CAB/SES is assumed 
to translate with the craft through its trajectory. The 
craft's axis system is also a right handed cartesian coor- 
dinate system with its origin at the center of gravity. 
„à The vector relationship between the craft's axes and the 


inertial reference frame are given by: 


X. cos y sin y 0 X 

А = sin y cos wy 0 Y 

Z 0 0 1 Ζ 
1 O | 


where the coordinate axis ín capital letters indicate vector 
quantities in the respective directions. 

Waves are also orientated with these two reference 
systems; the inertial frame to determine the wave properties 
and the craft's reference frame to determine the wave effect 
on the craft itself. Waves are assumed to be propagating 
in the direction defined by B, where B is the angle formed 
by the normal to the wave crests in the direction of propa- 
gation and the positive X axis. Waves propagate at an 


Enmplc relative. to the ship where y is defined by: 


Y=B- 4 2 
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For example, if the craft were not maneuvering, Y would 

equal B, since y equals zero (0), and for waves coming 

from the bow y would equal 180 degrees. | 
Considering a linear velocity Vn, a vector quantity 


E 
with components in the x, y, and z directions, 


ү. = 10 + ЈУ + kW (5) 
then, 
а μη EF, + κ ΣΕ, 


= m(iU + ДУ + kW) (6) 


Application of Newton's laws of motion with assumption of 
constant moments of inertía equal to rigid body values and 


small angle analysis leads to: 


Axial Force 
= - d 2 
КЕЛЕЕ at (mU) - mVR + mWQ (7) 
Lateral Force 
— = 
Y = ΣΡ. = TE (mV) + mUR mWP (8) 


Normal Force 


5 


2 = IF, = (mW) - mUQ + mVP (9) 


о, 


t 
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Rolling Moment 


аһ : š 
K = LM. шол = ТР (I,, ~ Ty OR - (R + PQ)I,, 
2 2 : 
+ (R" - Q mer + (PR - QI. (10) 
Pitching Moment 
dhy ° š 
M = LM, Rd > Ly Q N = п - (P + QR) Iy 
2 2 : | 
+ (Р -R Der + (PQ - R)I z (11) 


PUE. a SN Yawing Moment... E wr a: | E ο - 5 Е A ΜΝ .-:»... οσο ο DEL πα τον E» 
= dh, | | 
N = ¿N = w I,,R t (I. y - І) РӘ - (Q + RP)I 
2 2 к 
+ = = 
(Q P ту + (RQ Р)І,, (12) 
where 
η lxx -liy x P 
h SAET I zm | J 
y ух y y y2 Q (13) 
h, "lzx -Izy loz R 


and U, V, and W are linear velocities along and P, Q, 
and R are angular velocities about the x, y, and z axes 
respectively. Also associated with these equations are the 


kinematic relations 


фер + Usine (14) 
~ Q- ψοοςθ sinó 

8 = cos 8 (152 
y “kR Osing (16) 


^ cos80 cost 


20 





Next, consider the nonlinear inertial coupling as indicated 
in Equations (10) and (11). When the above equations were 

in the process of being formulated for the computer simu- 
lation, a computational approach was carried out using 
analog computers to determine the significance of nonlinear 
inertial coupling using the nonlinear model of inertial 
coupling terms. Results showed negligible difference in 
responses and the decision was made to simplify the equations 
of -τὋα by с ее ns Donc near ει coupling кете: 
— aic 5 Πο digital oee τ i 


Equations (7) through (12) may be simplified and written: 





U 0 0 0 0 0 X 
V m 0 0 0 0 Y - MUR 
W 0 m 0 0 0 Z 
= (17) 
В 0 0 Ix 0 "I. K 
0 0 0 0 M 
: y y 
| R 0 0 Pico 0 I N 


The above form of the state equations allows establishing 
the equations of motion in state variable form as a set of 
IESU order differential equations. W, P, and Q are 


augmented by the simplified relationships from equations 


W EZ | (18) 
Р = (19) 
Q= (20) 


2l 





. thus introducing additional state variables. These states 
may then be transposed back to the inertial axes by use of 
Equations (3), and trajectories are determined through 


simple integration procedures. 
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ТТ. PROBLEM STATEMENT 


The BELL 100B is assumed to be modelled by the computer 
program created by OCEANICS. The computer program is to be 
used to determine theoretically, the craft's turning 
characteristics in calm water and in simple waves. Turning 
characteristics relative to roll, pitch, yaw, and their 
rates of change are expected to compare favorably with 
available empirical data. 


. 
ow = 
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III. PROCEDURE AND DATA 


A. PHASE I DATA BASE 

Prior to development of the above problem statement and 
acquiring the necessary data, an initialized data base had 
to be established and trends, where trends existed, had to 
be noted. The data obtained in this study were obtained 
using the IBM 360/67 computer at the W. R. Church Computer 


Facility at the Naval Postgraduate School. 


i. Straight Runs 


= 494. 


* 
а a de ea? 
. 


a ae Bar ао о - 
longitudinal speeds (U) of 40 knots through 105 knots eis 
sive at 5 knot increments were used. Straight runs of 50 
seconds duration were made at each speed. Velocity was set 
at a constant value and thrust was allowed to vary as a 
function of surge only. For all speeds, except 95, 100, 
end105 knots, steady state values of thrust (T), pitch (0), 
and draft (d) were recorded and used as initialization values 
for a steady state run of 10 seconds. For the ten second 
run, all values remained constant; see Table II. 

Speeds of 95, 100, and 105 knots in actuality did 
not reach a steady state at the 50 second mark as did the 
other speeds. However, values of pitch, draft, and thrust 
were averaged over the last 20 second period. These values 
were in turn used in the 10 second run which again did not 
reach steady state. Thus, the values in Table II for these 


higher speeds are actually average values. Plenum pressure 
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Table Il. 


ORDERED SPEED (KTS) 


40 
45 
50 
55 
60 
"65 


70 


ως 


80 
85 
90 
95 
100 


105 


* 
* 


* 


*Average Values 


ONE ENGINE THRUST (LB) 


5091. 


Initialization Values 


5 


5490.2 


6002 


6589. 
7221. 
7872. 
8521. 
9107, 
| 9731. 
10254. 
10697. 
10839. 
11646. 


1276315 


. 


3 


25 


DRAFT (FT) 
17.212 
ПЕ 228 
1.166 
1.103 
1.003 
.956 


.874 


86 ᾿ 


.694 
. 598 
.498 
.362 
‚337 


° 322 


PITCH (DEG) 
.23 
. 26 
‚28 


.3l 


.39 
.39 
.40 
.57 
‚60 


„61 





was held at a constant 92.8 pounds per square inch for all 
runs. All succeeding simulation runs utilized these 
inltialized values, thus forming a steady state reference 
point for system perturbatlons. 
2. Straight Runs with Simple Waves 

To further the data basis for turning maneuvers, 1t 
was decided to make straight trial runs with varying wave- 
‘length (A), wave amplitude (a), and speed (U). Initial runs 
were conducted ш wave fronts сати from three different 
relative directions, 35 000° (astern waves), В OS (port 
beam waves), and B = 180% (bow waves). For each of these 
directions, speed (U) was varied (U = 40,50,60,70,85,90 and 
100 knots). Wave amplitude was set to a = 0.5 ft and wave- 
length A = 720 ft. Since this wavelength was ten craft 
lengths and the amplitude was relatively small, the situation 
approached that of an infinite wavelength. The initial 
conditions were set in accordance with data noted above. 

Table III shows the wave encounter period as 
experienced by ἘΠΕ Geralt ~the program computes initial 


encounter frequency in INCON as follows: 


= ος 
inter οπ(ρλ/οπ) (u cosB)/A . (21) 
Negative encounter periods occur when the craft overtakes 
the waves in a following sea. For waves from the bow and 
the stern, no motion was noted in the lateral plane. Pitch 


angle (8) did oscillate in both cases at the approximate 
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Table III. Initial Wave Encounter Period 


SPEED (knots) BOW WAVES (sec) STERN WAVES (sec) BEAM WAVES (sec) 


40 5.613 -105.26 πι ο 
50 4.960 - 30.342 11.859 
60 4.443 - 17.726 11.859 
70 4.024 - 12.520 11.859 
85 3.525 - 8,691 11.859 
90 3.385 - 8.887 | 11.859 
100 156 - 6.656 11.859 


e 





wave encounter period with slight phase difference. These 
trends deviated slightly at U = 100 knots. At 100 knots, 
still no motion was noted in the lateral plane. With the 
encounter period reduced to 2.136 seconds, venting, presence 
of air entry gaps along the sidewall, was noted at the 
center of gravity by the fact that draft was shown to be 
negative at various, though seemingly non-frequency related, 
intervals. The sidewalls of the craft did not completely 


Lege the water as noted by the c that the үш gauge 


. 
Πάει . 


. pressure did not fall below. zero as was seen ‚in later runs. 

Also noted in these runs where оет КЕ отап forces 
acted only in the vertical plane, was the fact that the 
relative difference between the accelerations of bow and the 
center of gravity never varied by more than 0.05 G in the 
"worst case! runs of U = 100 knots. In spite of this small 
variance, the craft's attitude in the vertical plane, 6, 
rose and fell with the rise and fall of the wave encounter 
period with both bow and stern waves. 

For л уез ап the port beam, the wave encounter 
frequency was the same encounter frequency as would be 
experienced if the craft were stationary. The bow and stern 
waves caused the initial velocity of U = 100 knots, 168.89 
ft/sec, to decrease to U = 161.69 ft/sec and U = 166.41 
ft/sec respectively by time t = 25 seconds. However, the 
beam waves produced only a slight oscillation about a 
nominal drop in speed of about .N ft/sec. Outputs for the 


beam wave case which vary significantly from the bow and 
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stern wave cases were pitch, roll, yaw rate, and displacement 
along the y-axis. It was interesting to note that an oscil- 
lation in pitch angle was produced when there were no 
explicit perturbation forces either in the z-direction or 
around the y-axis. The oscillation period appeared to be 
approximately that of the wave encounter period but slightly 
out of phase by one to two seconds. These oscillations were 
smaller than the oscillations noted in cases with bow and 


stern waves. 


, 1 Pes 


Further колты ion of data from the beam wave case 
showed that draft measurements oscillated at approximately 
the same frequency as wave encounter frequency but were 
180 degrees out of phase when considering maximum amplitudes.: 
Roll angle (4) and yaw rate (R) also showed a marked perio- 
Gicity with respect to the encounter frequency. However, as 
speed (U) was increased this relationship to the encounter 
frequency seemed to diminish. Roll angle showed no discern- 
ible frequency whatsoever at 100 knots. | 
The next Buen taken was to reduce the wavelength of 
the simple wave. In the same manner as explained above, 
with seas om three directions, computer runs were made 
using a wavelength A = 360 feet or five craft lengths, with 
speeds U = 40, 50, 60, 70, 85, 90, and 100 knots. Results 
with the following exceptions were approximately the same 
as the above case: even at the 100 knot level the pitch 


angle showed a definite relationship to the wave encounter 


frequency oscillations; venting did occur for the bow waves 
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and the stern waves but the excursion in plenum gauge 
pressure was much greater in both cases; differences in 
bow and center of mass accelerations were greater; longi- 
tudinal velocity dropped off more in both the bow and stern 
seas case, but longitudinal velocity for the port beam seas 
show a slight but almost constant decrease, vice an oscil- 
lation over the 25 second problem period; maximum pitch 
excursions for the beam seas were about the same in both 
cases, +0.5 degrees, but this ur on was ο ο 
_ less for the bow and stern sea cases. | | 2 
The next step іп setting En basis was ИЯ 
parameter variances over a wide range of representative 
values. To accomplish this, it was decided to utilize 
wavelength values of relative craft lengths. Thus, wave- 
length (X) was varied as one-half a craft length, A = 36 
feet, one craft length, A = 72 feet, two craft lengths, 
A = 144 feet, and three craft lengths, A = 216 feet. Іп 
that the plenum average height is slightly over six feet 
and that drafts of 1.5 feet have been measured with only 
a 0.5 foot wave amplitude, wave amplitudes of a = 0.5, 1.0, 
and 2.0 et were chosen as representative values.  Further- 
more, it was noted that in the program, the input wave 
amplitude "a" was in fact an amplitude, but when considering 
effects of sea waves, the peak to peak value, the crest to 
trough value or the wave height is the measurement usually 
considered. In other words, the above values of wave 


amplitudes chosen have heights of 1, 2, and 4 feet 
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respectively. The latter height leaves little room for 
heave mctions and therefore, when water contacted the top 
of the plenum, whenever fan horsepower went negative, or 
whenever the plenum gauge pressure went negative, the run 
was considered to be aborted. 

The same procedure, as followed for the above cases, 
where A = 720 feet and A = 360 feet, was utilized for the 
varying wavelength and wave amplitudes. For each value of 

wave une a = 0. le 0, and 2. gn a total of 84 simu- 
lations were run ‚ varying wavelength, wave „direction, and 
speed. An example of the data collected is shown for beam 
waves, in Table IV. Data for bow and stern waves were also 
collected, but the data were omitted due to later usage of 
beam waves as an initial condition for turns. 

Behavior of the CAB/SES in the above wave conditions 
was best compared by categorizing the waves by wave direction 
and not by wave amplitude. The most "well behaved" of all 
the wave directions was the beam wave. Well behaved in this 
case.means that of all the wave directions, beam waves 
provided fewer run abortions than bow and stern waves. Of 
the 36 Parise Окс δα in Table IV, only three runs were 
aborted due to water contact with the plenum deck. The bow 
and stern waves had Bas aborted runs each. 

For bow waves pitch angle maximum excursion was 
greatest for the longer wavelength, lower speeds, and higher 
wave amplitudes. The maximum pitch for bow waves was 


recorded at t = 20.5 seconds of the run using the wavelength 
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А = 216 feet and amplitude a = 2 feet at speed U = HO knots. 
The frequency of oscillation of the pitch angle seemed to 
correspond roughly to the wave encounter frequency. However, 
in that print output was a one-half second interval, 1t was 
difficult ot determine phase delay times, and for shorter 
encounter periods it was impossible to determine any phase 
or frequency relationship at all. Starting with the one 
foot amplitude waves and continuing through, no pattern 
could be determined with gauge pressure fluctuation for the 
bow waves with the Е πώ να that as speeds 
IN ceased, makinan I DAY eier Я This 
was also the case for stern waves. Maximum acceleration, 
pitch producing, differentials appeared between the center 
of gravity and the bow increasing with wavelength until 
À = 144 feet and then dropped off considerably for smaller 
wavelengths. Interestingly final speeds of runs which were 
not aborted hit a low point at A = 144 feet. 

Stern waves, aside from producing the same number 
of aborted runs as bow waves, produced decidedly more 
adverse reactions ie craft. For all wave amplitude 
groupings, A = 144 seemed to affect the craft more adversely 
than any other wavelength. Maximum pitch and pitch excurslon 
were produced with wave amplitudes of a = 0.5 and 1.0 feet 
and speeds of 40 knots. Additionally, the most violent 
plenum pressure gauge fluctuations, bow and center of mass 


accelerations differences, and speed loss were also noted 
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for А = 144 feet. As one would expect, there were no yawing 
or roll moments from either the bow or stern waves. Aborted 
runs due to negative fan horsepower seemed to occur whenever 
gauge pressures went above 120 pounds per square foot. 
Results for the beam waves are presented in Table IV. 
Generally ος. maximum pitch and pitch excursions 
increased with speed amplitude and decreasing wavelength with 
only' two exceptions for non-aborted runs during the a = 2 


foot simulations. The frequency of oscillation seemed to 


. _. Match wave encounter frequency in, every case with little, or 


in some cases, seemingly no phase difference. Acceleration 
variances between the bow and the center of mass were 
considerably less than values noted for either bow or stern 
wave cases, as were plenum gauge pressure fluctuation 
extremes. The low point for final speed loss in most cases 
appeared at values of A = 72 feet. The major difference 
between bow and stern waves E beam waves, for the non- 

' maneuvering case, were the lateral component wave forces 
which produce roll and yaw moments. For the beam wave 
simulations, roll in every case decreased with speed. Also, 
values ОС o1] increased up to A = 72 feet and then dropped 
off at À = 36 feet for cases where speed was held constant. 
Both yaw rate and roll frequencies of oscillation were 
related to encounter frequency with varying phase time lag. 
The yawing effect was noted to be considerably higher for 
higher amplitude waves and then dropped off within amplitude 


groupings when moving from A = 72 feet to A = 36 feet. 
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Bee PHASE TI TURNING DATA 
1. Constant Flat Turns 

To provide a wide basis for turning data speeds of 
U= 40, 50, 60, 70, 85, 90 and 100 knots were each simulated 
Pomoc vurms Of constant 5, 7.5. 10, 12.5, 15, 17.5, and 
20 degree rudder. For each rudder used under each speed the 
problem was allowed to run on a straight course for two 
seconds prior to moving e ruder. Initial condition values 
s determined above were not varied. At the two second mark,. 
the rudder was moved” five degrees’ per second until the final 
rudder was reached. Thus, the final rudder angles were 
reached at times t = 3.0, 3.5, 4.0, 4.5, 5.0, 5.5, and 6.0 
seconds respectively for the above rudder angles. Figures 
1 through 9 show roll, pitch, and yaw for representative 
speeds U = 40, 85, 505 100 knots. In all cases the ordinate 
units were measured in degrees and the abscissa's units 
denoted problem time up to 25 seconds. To reduce confusion 
on the plots, only 5, 15, and 20 degree rudder values were 
plotted. For these turns it was noted that values of pitch 
and roll tended to oscillate about an average value. As 
time increased, the average value of roll/heel tended to 
decrease. This may be seen in Figures 1, 4, and 7, where 
the negative values of roll indicate a heeling motion to 
port for 5, 15, and 20 degree turns at 40, 85, and 100 knots. 
In the analysis of the yaw rate plot, Figure 10, yaw rate 


also oscillated around an average negative value, indicating 
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a turn to port. Yaw rate appears to be constant when 
considering the yaw versus time curves in Figures 3, 6, and 
9. 

In all cases, there was a steady decrease in speed 
throughout the 25 second turn problem; Table V. The more 
rudder used, the more final speed was down from the initial 
value. Average final roll angle, or heel, increased with 
“increasing speed; Table VI. Yaw and yaw rate reacted in 
like manner but also increased with rudder angle as ' 
566 bles Vit and VIII. Pitch angle Was slightly more, 
erratic. For 40 and 50 knot runs, final average pitch 
angle for five degree turns were slightly less than the 
initial condition steady state pitch angles. In both cases 
as rudder angle was increased, pitch angle decreased. For 
all other speeds, the five degree final average pitch angle 
was higher than the initial condition value for the corres- 
ponding speed. As larger rudder angles were used the final 
average pitch angles increased and then dropped off at 
rudder angles of 7.5 degrees for 60 knots and 85 knots, and 
dropped off at 10.0 degrees rudder for 7/0 knots and 100 
knots; Table IX. Turning radii varied considerably for 
five degree rudders on a speed to speed basis. However, | 
turning radii seemed not to be dependent upon speed for 15 
degree and 20 degree rudders; Table X. 

2. The Roll Phenomenum of the HO Knot Turns ` 
Whiter plovtine various turn data for the forty knot 


turns versus problem time, it was noted that values for roll 
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RUDDER ANGLE IN DEGREES 


15 


20 


10. 


12% 


15. 


I7. 


20. 


Table V. 
40 50 
66.7 83.0 
60.7 74.4 
56.4 68.7 

Table VI. 
К: 

352.593 

.469 .766 

. 541 .872 

.577 9897 

. 580 .871 

.547 .824 

.496 «757 


Final Time Simulation Speed (feet/second) 


SPEED IN KNOTS 


60 
99.1 
87.1 


ПЭ 


Average Final 
SPEED IN KNOTS 
. 60... 


.86 


1.10 


1.24 


1,25 


1.20 


REST 


297 


70 
114.7 
98.5 


89.6 


το 
1.17 
1.49 
1.62 
1.62 
1.54 
1.42 


1023 


5; 


85 90 100 
136.9 143.1 153.2 
112.5 116.4 123.5 


101.0 104.2 109.8 


Roll Angle (degrees) 


85 90 100. . 


TS 


12955877205 772,50 
2.192: 372.33 72,50 
2.17 2.34 2:63 
2.03. 2.19 2.47 
12932271796 557520 


1.60 1.71 1.90 





RUDDER ANGLE IN DEGREES 


15 


20 


10. 


12. 


15. 


17. 


20. 


Table VII. 


40 
16.4 
53.7 


70.6 


40. 


1.0 
1.4 


1.8 


SPEED IN KNOTS 


50 60 
20.5 23.3 
64.8 1533 
85.3 

Table VIII. 


70 


26.3 


84.6 


80 


28.9 


Final Yaw Angle (degrees) 


90 


29.9 


100 


33.0 


97.8 101.0 109.0 


99.0 111.4 ‘129.6 134.5 145.5 


SPEED IN KNOTS 


DOF 2. 


М 


‚8 


2.1 


sm 


3.5 


60 


1.0 


ς > ° 


PESE 


40 


70 


Average Yaw Rate 


18 


2.6 


33 


5105 


5.1 


(degrees/second) 


100 i 


% 

1° 13 
1.8 2.0 
2.6 2.8 
3.3 3:9 
4.0 4.3 
4.6 4,9 
5.2 ЭЭ 





RUDDER ANGLE IN DEGREES 


os, 


10.0 


1,225 


15.0 


17.5 


20.0 


Initial 


Il» 


20 


40 

.220 
.210 
πο 
. 178 
.154 


125 


..095 


„25 


40 


5800 


1590 


1130 


Table IX. 


50 


„279 


‚276 


‚267 


‚250 


‚218 


“1/1 


К 


‚28 


Table 


50 


5850 


1600 


1120 


X. 


SPEED IN KNOTS 


60 


. 340 


‚351 


. 350 


.328 


.298 


„252 


.206 — 
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SPEED IN KNOTS 


60 
6000 
1620 


1140 


41 


70 

.407 
.440 
.445 
.428 
.387 
.326 
0253 


‚36 


70 
6100 
1640 


1130 


85 


‚603 


‚652 


„590 


„279 


‚526 


.454 


.363 


‚39 


Average Turn Radii 


85 


6750 


1640 


1130 


Average Pitch Angle (degrees) 


90 100 
.646 «673 
.728 | .740 
.706 ПОТ 
‚626 „747 
.560 «658 
.489  .568 
.393. ο... 
.40 62 
(feet) 

90 100 
6720 6750 
1650 1660 
1140 1135 





Ἱ πο που ο ο  υυτς wichi expected results; see Figure 1. 

The value of roll for the first positive cycle, roll to 
starboard, for five degrees was greater than the value of 
the twenty degree turn by about 0.8 degrees. Furthermore, 
throughout the rest of the simulation, the amplitude of 

the five degree rudder oscillation was substantially greater 
than the amplitude for the twenty degree rudder. | 

To gather data to explain this apparent phenomenon, 
initial results were first double checked both with smaller 
tolerancés for the: Runge-Kutta-Merson: numerical. integration 
technique and a constant step size; see computer program 
section. Results showed initial calculations to be correct. 
Next, values of roll were retabulated for rudder angles of 
five, fifteen and twenty degrees; Table XI. It was noted 
that the values for the fifteen degree rudder were the same 
point by point as the five degree rudder until time t = 3.5 
seconds; the rudder had stopped moving for the five degree 
Е вме Т = Sececomds and time 3.5 Seconds the rudder 
was still swinging for the fifteen degree rudder simulation. 
Likewise, the values of roll for the twenty degree turn 
were the same at time t = 5.5 seconds as the values for 
the fifteen degree turn. 

The program was then modified, see the computer 
program section, to provide an optional output which included 
summations of axial ON lateral forces (Y), vertical 
forces (Z), roll moments (K), pitch moments (M), and yaw 


moments (N); see Figures 11 through 16. Further modification 
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was made to the basic program to provide the components of 
the roll moment (K). These include the bowseal, the stern- 
seal, the rudder, the side walls and appendages, the pro- 
pellers, the bubble, as well as components from wave and 
aerodynamic inputs. As above, with roll angles, the summa- 
tion of moments about the x-axis were equal through time 
3.0 seconds. At t = 3.5 seconds, the summation of noments 
for the five degree turn case jumped to over 44,000 ft. lb. 
where as the summation of moments was only about 14,000 ft. 
- 1b: for the twenty degree turn. See Table XII. for precise 
tabulation results of times t = 2 seconds through t = 7 
seconds for the summation of moments about the x-axis. 
ο Turns with Simple Waves 

The final steps of this procedure included deter- 
mining what and how many simulation runs to make to insure 
valid data. Original intentions were to establish a field 
of various wave lengths, speeds, and wave amplitudes from 
all three wave directions as utilized in the above procedure. 
However, during those runs it was noted that the craft's 
reactions to parametric changes were highly dependent upon 
wave encounter frequency. The craft reacted in a much 
more seaworthy manner with beam waves. For higher veloci- 
ties and amplitudes and smaller wavelengths, simulation 
runs were often aborted due to loss of plenum pressure, 
water contact with the deck of the plenum or due to negative 
plenum fan horsepower caused by high pressures within the 


plenum. It was noted that the craft will experience 
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Table XII. Summation of Principle Roll Producing Moments 


Five Degree Rudder Simulation 


TIME M M M M RUDDER 
x rudder bubble sidewall 
(sec) CEEC Ib.) СЕВЕ) (ft .1b.) Ga Tb.) (deg) 
2.0 0.0 0.0 0.0 0.0 0.0 
2.5 32252250: 82. 011.0 “ 1700.0 10,359. 2 
3.0 19,994.0 63,742. 4872. 48,782. 5.0 
345 44,603.0 36,589. 20132; 91,491. 5.0 
4.0 49,689.0 28,959. 20510 89,133. 5.0 
4.5 10,847.0 29,719. 9673. 45,510. 5.0 
5.0 29,839.0 32,491. 3517. 854. 5.0 
D» 45 134.0 35.1273 4501. 14,864. 5.0 
6.0 28,444.0 37,465. 4676. 6,630. 5.0 
6.5 5,060.0 39,503. 7010. 48,214. 5.0 
2.0 32,804.0 40,868. 17522 T 6 2; 627. 5.0 
Fifteen Degree Rudder Simulation 
2.0 SAME AS ABOVE 0.0 
2.5 255 
3.0 5.0 
3.5 14,672. 78,945. 21,788. 104,762 7.5 
4,0 37.159. 95,386. 29,625. 147,688 10.0 
4.5 28,570. 119,132, 28,287. 158,146 12.5 
5.0 12988, 138,015. 19,880. 146,565. 15.0 
5:5 Е 120,520. 10,047. 1220223, 15.0 
6.0 21,690. 119,312. 1,469. 98,770 15.0 
6.5 30,961. 121,867. 1,459. 90,146. 15.0 
7.0 18,247, 126,382. 3,159. 107,645. 15.0 


_ Twenty Degree Rudder Simulation 


20 SAME AS ABOVE 0.0 
2.5 209 
370 5.0 
3.5 y 
4.0 10.0 
4.5 12.5 
5.0 1,580. 102.970. 15.0 
5.5 28,056 187,987. 11,593 140,705. 17.5 
6.0 31,414. 1705725. 9,062. 101,486 20.0 
6.5 120922, 168,321. πο το. 190,036. 20.0 
1.0 21,524. 170,424. 14,910. 197,390. 20.0 


NOTE: All moment producing elements are not included. 
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varying encounter frequencies with any initial wave direction 
utilized since wave encounter frequency is a function of 
wavelength, longitudinal velocity, and relative wave 
direction. For each simulation where the ship is maneuver- 
ing an instantaneous encounter frequency is determined. 
Thus, it was decided to start each problem with waves from 
the port beam, allow the craft to follow a straight course 
for the first two seconds and then put the rudder over at 
time t = 2 seconds with a rudder swing of five degrees per 
second as was done above for the flat turns. Again initial- 
а о т пе ео оса овса 
шеге used. | | o | L | 
At this point another cecision had to be made 1п 
order to provide enough useable data. The original plan was 
to make sets of twelve computer runs for each of the speeds, 
О = 40, 50, 60, 70, 85, 90, and 100 knots. Six-hundred- 
seventy two simulations averaging seven minutes each were 
mearly out Of the question. On a second cut, it was 
decided to use each of the above speeds with eight runs. 
Bach run consisted of seven cases with varying rudder angles 
of five degrees through twenty degrees in 2.5 degrees 
increments. The eight runs used various combinations of 
previously used wave amplitudes a = 0.5, 1.0, and 2.0 feet 
and wavelengths A = 216, 144, 72, and 36 feet. For each 
group of runs, A = 144, a = 2; AX = 72, a= 2; 4 = 144, a = 1j 
and À = 72, a = 1, was deleted. These values were deleted 


not only because of the computer time requirement, but also 
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due to the fact that these four sets provided repetitive 
values of wave amplitude to wavelength ratios: 


Max slope = 360° · t (22) 


as determined in the subroutine INCON. 

Table XIII is presented to show the format for data 
collected at all speeds. Particular turning maxima and 
minima are noted in Table XIV. The terms Ad and AB in 
Table. XIV are maximum roll and pitch excursions respectively. 
‚Maximum values were absolute values while the minimum values 
were the minimum of the maximum values of collected data. 

In general, roll seemed to increase with decreasing wave- 
length until À = 36 feet. At this point roll values 

decreased significantly. In all cases, roll values were 

less for wavelengths of 36 feet than they were for wavelengths 
of 216 feet. 

Pitch generally decreased from A = 216 feet with 
decreasing s. stis: Side slip angle, the angle formed 
by the direction of force summation and the ship's longitu- 
dinal axis, varied with increasing rudder angle. Rudder 
angles of five degrees, fifteen degrees, and twenty degrees 
corresponded to average side slip angles of 2.3, 5.8 and 
7.3 degrees respectively. Yaw angle and yaw rate increased 
with increasing rudder angle and speed. Further observations 
not directly associated with turning data were gauge pressure 


readings and final speed relationships. Again, as with 
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Table XIV. Maxima/Minima Summary for Turns with Waves 


MAXIMA MINIMA 
Value Simulation Value Simulation 
(deg) U a А ги (deg) U a A ὃ τιιά 
(ft/sec) (ft) (ft) (deg) (ft/sec) (ft) (ft) (deg) 
ROLL 4.63 40.0 0.5 72.0 5.0 1.68 40.0 0.5- 36.0 5.0 
ROLL 9.35 40.0 0.5 72.0 5.0 0.90 100.0 0.5 36.0 5.0 
PITCH 3.09 40.0 1.0 216.0 5.0 0.66 50.0 0.5 36.0 5.0 


PITCH 4.51 50.0 1.0 216.0 15.0 0.27 100.0 0.5 36.0 15.0 
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straight runs with waves, runs were aborted when plenum 
gauge pressure surpassed 120 psf. Runs with wavelengths of 
36 feet provided a slight decrease in gauge pressure. Speed 
drop-off was also noted with increasing time, speed, wave 
amplitudes and decreasing wavelengths with the exception of 


36 foot wavelengths where final speed values increased. 
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ТУ. DISCUSSION/EVALUATION OF DATA AND PROCEDURE 


The original plan of attack for obtaining turning data 
was both broad and sound considering the fact that no 
information other than Reference 1, the user's manual, was 
available for simulation information. The data base was 
established by determining relationships for straight runs 
Bach simple waves, and by determining craft reaction to 
Gaim water condition turns. The procedural part of the 
original plan was to utilize a wide set of dote points by, 
varying speeds from 40 knots to 105 knots, varying. rudder 
angles from five degrees to twenty degrees, varying wave- 
lengths from one half of a craft length to ten craft lengths, 
varying wave amplitudes from 0.5 to five feet, and varying 
wave directions from 180, 90, and 000 degrees relative. 

The sheer magnitude of possible data accumulation demanded 
that data points De minimized and chosen so as to provide 
Wide spectrum of relevant data. Thus, after initial 
conditions had en found for speeds in five knot increments 
from 40 knots to 105 knots, the field was cut to only seven 
Speeds for the straight runs with waves. Prior to proceeding 
to the turns with waves, it was decided to limit initial wave 
direction to port beam waves. The beam wave data from the | 
flat turns were further subdivided and only 40, 60, and 

100 knot information was extracted. Finally, of 84 possible 
Simulation runs for each of the seven speeds for the turns 


with waves simulation, 28 of these were deleted and only a 
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total of 24 runs, 5, 15, and 20 degree rudder turns, were 
recorded for the broad analysis. In the final analysis, 


only the 40 knot turns were presented in this study. 


A. PHASE I DATA BASE 
1. Straight Runs 
The straight simulation runs with no waves set the 
initial conditions for all succeeding runs. Although steady 
` state conditions were not determined for the higher speeds, 
U= 95, 100, and 105 knots, it was decided to accept the 
average values as steady state conditions. 
‘2. Straight Runs with Waves 
The first simulations with waves utilized long 
wavelengths and small wave amplitudes with the intent of 
extracting information due to nominal perturbations of the 
wave effects at what approximated a wave of infinite length 
relative to the craft's dimensions. Unfortunately, the 
first wavelength used, A = 720 feet, yielded more reaction 
than had been expected. Rather than Ее wavelength 
to approach the originally desired relative infinite wave- 
ЖО орЫН. it was decided to halve the wavelength to five craft 


lengths or i = 360 feet. The final decision to use wave- 


lengths of A = 216, 144, 72, and 36 feet was based on two 
main ideas. Firstly, there had to be a choice of wavelengths 
that would provide relevant data. Secondly, actual seas 
encountered by such a craft would not usually be of the 


extreme wavelengths as encountered in open seas by larger 


vessels. 
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matche procedure section above, it was μη that 
Wave encounter frequency played an important role in the 
ship's stability characteristics. Roll, pitch, and yaw 
all oscillated with some relationship to this wave encounter 
frequency. The relatively gross one-half second print inter- 
val used to determine the values for the beam wave simula- 
tions showed only a general relationship and negligible 
phase information. Figure 17 shows wave amplitude for a 
straight ДО knot run with A = 216 and a = 0.5 feet. When 


wave amplitude was compared with Figures 18 and 19, plots 


". .Of roll and pitch, it was. readily seen that. roll and pitch 


were definitely wave encounter frequency related with only 
a Slight delay time. Interestingly, the roll and pitch 
producing moment plots, Figures 20 and 21, showed no such 
clarity for frequency or phase relationships due to the 
irregular formation of the moment curves. However, when: 
the roll and pitch curves were compared to their respective 
moment curves, it was noted that they were approximately 
180 degrees or about one-half of a wave encounter period 
out of phase. 

The loss of speed noted in the wave procedure was 
due to the increased drag caused by the wave. Figure 2? 
shows wave drag in the x-direction with varying rudder 
applications under calm water conditions. The important 
poet of higure 22 for this example is that it illustrated 
an initial wave force which was constant until the rudder 


was put in motion at the two second time of simulation. 
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It is this initial wave drag which when coupled with the 
additional hydrostatic drag forces caused loss of speed. 
As velocity decreased more wave drag was experienced, see 
reference [2], and speed was further decreased. The end 
result was that eventually a steady state would be reached 
whereby the initial thrust as set into the simulation for 
the desired initial speed would eventually equal the total 
drag forces thereby reaching a steady state speed. Times 
Tor ο steady state и were У longer 
than the 25 second simulation runs used in this report. 

For the bow and stern waves the roughest perturba- 
tions experienced by the craft were simulations with А = 144 
feet. Beam waves similarly provided maximum distrubances 
at А = T2 feet. The cause of this had to be attributed to 
the craft's dimensions. For longer wave lengths the craft 
experienced a rise and fall with the wave encounter fre- 
quency. However, for wavelengths of equal size to or smaller 
than the dimension in question, the 75 foot length for the 
bow and stern waves and the 36 foot width for beam waves, 
angular displacement excursions were greatly reduced due 
to the rapidly changing wave forces. The fact that these 
angular excursions were less was not in itself a fortuitous 
Peeurrance. As noted in Table IV, the only aborted runs 
were at the lower, 36 foot, wavelengths for the beam waves. 
EMUuTtlonally, all of the aborted runs for bow and stern 
waves, not shown, occurred at wavelengths of 72 feet or 


less. The reason for these aborted runs was the excessive 
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angular accelerations experienced by the craft due to 
increased wave encounter frequencies and increased wave 


amplitudes. 


B. PHASE 11 TURNING DATA 

1. Calm Condition Turns 

Most of the data collected did not indicate any 

unexpected results. Roll, pitch, and yaw rate all exper- 
jenced some type of decaying oscillation commencing with 
Козге С = 2 seconds, the beginning of rudder motion.” This 
[udder motion was considered © be а ramp ao ипе. Uus 
meached its final desired С and Ее a ent 
at that value throughout the simulation. With the cessation 
n rudder motion, force values due to rudder change also 
began an oscillatory decay. Roll characteristics will be 
discussed in greater detail below, but for the most part 
roll attempted to reach a steady state heeling attitude to 
port. Pitch final values, as noted in Table IX and illus- 
trated in Figures 2, 5, and 8, for 5, 15, and 20 degree 
turns with speeds of 40, 85, 100 knots, oscillated about 
values less than the initial condition values. Figures 2, 
5, and 8 indicate that the more rudder’ applied the greater 
the average pitch angle. Yaw as opposed to roll and pitch 
seemed to be fairly steady when viewing Figures 3, 6, and 
9. However, Figure 10 shows that a 40 knot simulation did 
experience oscillations in yaw rate producing more oscillation 


with the smaller rudder angles. 
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2. The Roll Phenomenom of the 40 Knot Turns 
As noted in the procedure, positive peak overshoot, 

roll to starboard, was much greater for five degree rudder 
than for either 15 or 20 degree rudders; Figure 1. The 
whole problem was eventually resolved by calculating the 
breakdown of the roll producing or the components.of the 
moments about the x-axis. When data such as found in Tables 
x and XIT were properly correlated on a time basis, it 

was seen that ш У the а. were exactly е. 
until КО (осел, er the er earth its ο. И ο 
mu Limes greater. than three seconds for the five degree 
turn, and times greater than five seconds for the 15 degree 
turn, a different problem was generated than for the 15 
degree and 20 degree turns respectively. In that the clock- 
wise motion of the rudder produced greater negative roll 
moments with increasing rudder at print time 3.5 seconds, 
the negative moment component due to the rudder was much 
greater for the 15 and 20 degree rudders than for the five 
degree rudder. ит. at print time 5.5 seconds, the 
negative moment for the 20 degree rudder was significantly 
larger than the moment for the 15 degree rudder. As the 
craft begins to turn the dynamics are such that the summa- 
pron of forces, and for that matter the direction of motion, 
15 not always directed along the craft's longitudinal axis. 
jm. fact, for a turning maneuver with constant rudder a new 
ime of force is generated. This direction of motion 


relative to the craft's heading is known as side slip angle. 


58 





In the simulation program, side slip angle is simply 


determined by 


els 


Side Slip Angle = (23) 


in radians. The initial motion of the rudder provides this 


| panes E Ed CC LODS At ы Point, zen) ne Со: 


it was E s ee ΠΠ A Кс апше ne cc 


„component of roll moment due to the sidewalls increased 
until the moment component peaked within πα of a 
secónd after the rudder motion had stopped: Average side- 
slip angles were 2.1, 5.8, and 7.2 degrees for 5, 15, and 
20 degree rudder turns respectively. Obviously, the in- 
creasing side slip angle produced increased moments with 
higher rudders, however; as noted above, these values did 
not peak until after cessation of rudder motion. Thus, at 
time 3.5 seconds the slightly large positive moment due to 
the sidewalls coupled with the notably smaller negative 
moments for the rudder produced a much more positive roll 
reaction, roll to starboard, than in the cases of the 15 
and 20 degree turns. 
3. Turns with Simple Waves 

ος πο πε ofthe first runs of turas with waves 
showed that the best conditions generated by the beam wave 
case soon deteriorated as the craft's heading turned toward 
the wave front. Due to the one-half second print interval 


decided upon for all the runs, frequency and phase analysis 
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charca cLeristies were not carried out. Figure 
23 plots the wave amplitude for a 40 knot 15 degree rudder 
run with A = 216 ft and a = 0.5 ft. Notice in particular 
how encounter frequency had increased near the end of the 25 
second run. Figures 24 through 26 show the same 40 knot 
turning run compared to a 40 knot turn without waves in 

roll ΕΡΕ 22s xn IE env. Notice that the frequencies 
jor the doting. о ποπ... E Some dégree the": 
ο ο ος ΞΕ ao no waves with only slight 
time delay. Notice "IS that in Figure 26, the turning a 
. path for the wave.case seems to have a much smaller radius 
mnan the radius for the no wave case. This is attributed 
to the z-axis moment added by the wave. Although cases 
were not run past the 25 second mark, it is assumed that 
the turning radius for the wave case would be correspondingly 
larger when the z-axis contributed by the wave acts against 
the direction of turn when completing a 360 degree turn. 

As with straight run waves, turning runs with waves 
seemed very ορ. on wavelength. Slower runs with long 
wavelengths yielded maximum perturbations in roll, pitch, 
anar roll rate. For turns with waves, the drop off point 
for some parameter excursions, and the minimum for others 
was for wavelengths of A = 36 ft. This seems reasonable 
since as in the above straight run simulations with beam 
waves, the beam dimension of 36 ft produced the same results 


when considering wavelengths. 
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V. CONCLUSIONS 


A. TURNING CHARACTERISTICS 

The study of any craft's turning characteristics in 
the final analysis must be based upon actual turning data 
generated by a prototype craft. For the purposes of this 
as correlating data" éould:háve" generated: more ‘definite +, 
conclusions as to the validity of turning characteristics 
presented above. ee since no such data has been made 
available only generalizations can be made as to the 
Simulation capability of the Oceanic's en 

Basically, turning characteristics seemed to be compatible 
with previous knowledge and past experience. From a controls 
point of view, data tends to indicate that with this 
particular CAB/SES as modelled would be optimally utilized 
at higher speeds using larger rudder angles. Of course, the 
dangers of high speed maneuvering are obvious when considering 
increased and varying wave encounter frequencies. In that 
the craft was ЕОР for high speed operations, control 
efforts should be directed at reducing both the acceleration 
and the maximum excursions of attitude parameters. As it 
is, the simulation showed only such small changes in these 
parameters that if the same perturbations were experienced 
by a displacement vessel, they would in general go unnoticed. 
However, since the CAB/SES is not a displacement vessel, 


these stability parameters are much more important. 
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ral the simulation seems to be accurate when 
considering horizontal plane trajectories. As far as the 
above simulation runs were concerned, data obtained for 
vertical plane motion, heave in particular, leave room for 
some doubt. Original runs used the variable Runge-Kutta- 
Mersen numerical integration. As ade were introduced 


to sr ο ο zuns; more ut en и was u due 


. . a ο 
ο at pus ` ` < = i. 


to s small ratio ЕО et ЕЕ d E numerous 
step size doubling and halving iterations. To prevent 


excessive ld Di the tolerances in some cases were 


| ea осу Е басне small error. . Most. trajectories 


were relatively unchanged; however, the most sensitive 
output parameters dealt with plenum pressures and motions 

in the vertical uen Final runs with turns and waves 

used a fixed step size of 0.005 seconds which was determined 
to be the average step size used by the Runge-Kutta-Mersen 
method. At the time of this writing, it seems that more 
accurate results could have been obtained with a 0.001 
second step size while sacrificing four times as much 
computer time. For the purposes of this study, it is felt 


that the .005 second step size was sufficient. 


B. RECOMMENDATIONS FOR FUTURE STUDY 

The first and foremost objective is to determine if the 
loads and motions simulation is valid., When this study was 
initially undertaken it was with the idea of eventually 


comparing simulation data with actual test data. At the 
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time of this writing test data had not been made available. 
Once these data are compared more confidence will be gained 
in the loads and motions simulation thus enabling modifica- 
ЕЛЕ ипе simulacion Tor larger craft, Additionally, 
steps should be taken to modify the simulation for the XR-3 
to further develop understanding and confidence in the pro- 
gram. Once these νους have been taken the u obvious 
a is to model lee ου στα ton gross E 'CAB/SES. 
Once it has been proved that the basic program is 
adequate ano ο. IC car be modified, a very interesting 
Д field Or Bed would be an optimal Control system: i Items 
ENrtroCcondrepcusror These Controls aside fron standard 
rudder and propulsion apparatus would be fan pressure 
control for the seals as well as the plenum, fins, and any 
other type of control surfaces that would enhance the 
stability control that would be required for high speed 


operations. 
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Figure 2. 


Forty Knot Turns: Pitch Angle Versus Time 
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Yaw Angle Versus Time 
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Eighty-five Knot Turns: 
Roll Angle Versus Time 
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Eighty-five Knot Turns: 
Pitch Angle Versus Time 
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Figure 6. Eighty-five Knot Turns 
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One Hundred Knot Turns: 
Roll Angle Versus Time 
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Figure 8. One Hundred Knot .Turns: 
Pitch Angle Versus Time 
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Figure 9. One Hundred Knot Turns 
Yaw Angle Versus Time 
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Figure 10. Forty Knot Turns: 


Turn Rate Versus Time 
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Forty Knot Turns: 
Summation of Axial Forces Versus Time 
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Figure 12. 


Forty Knot Turns: 
Summation of Lateral Forces Versus Time 
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Figure 13. Forty Knot Turns: - 
Summation of Vertical Forces Versus Time 
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Forty Knot Turns: 


Summation of Roll Producing Moments 
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Forty Knot Straight Run with Waves: 
Wave Amplitude Versus Time 
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